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A quantum field theory approach is put forward to generalize the concept of classical spatial light 
beams carrying orbital angular momentum to the single-photon level. This quantization framework 
is carried out both in the paraxial and nonparaxial regimes. Upon extension to the optical phase 
space, closed-form expressions are found for a photon Wigner representation describing transforma- 
tions on the orbital Poincare sphere of unitarily related families of paraxial spatial modes. 
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I. INTRODUCTION 

Photons are among the main carriers of information. 
This information can be encoded in their energy, linear 
momentum, and polarization state. In recent years, an- 
other degree of freedom for photons has been recognized: 
their orbital angular momentum (OAM). In 1992 Allen 
and coworkers |1| showed that optical paraxial cylindri- 
cal beams having an azimuthal phase dependence of the 
form exp(il(f)) carry a discrete OAM of lh units per pho- 
ton along their propagation direction. This angular mo- 
mentum produces a mechanical effect (induces a torque) 
when suitable light patterns interact with matter; it can 
be transferred from spatial beams containing phase dis- 
locations on their axis (e.g. optical vortices) to suitable 
trapped microscopic particles in optical tweezers [2J. 

At the quantum level, considerable interest has been 
brought for quantum information processing exploiting 
single and entangled photons prepared in a superposi- 
tion of states bearing a well defined OAM [1 H U H, H 
H, 0, 0. Indeed, one of their main distinguishing fea- 
tures is that, at variance with two-level quantum states, 
or qubits, OAM photon eigenstates involve the more gen- 
eral case of d-level (d > 2) quantum states or qudits. 
Such a generalization to multidimensional states allows 
to extend quantum coding alphabets without the need 
to increase the number of entangled photons, providing 
also a more secure quantum cryptography. Since fewer 
photons are needed, the multidimensional approach re- 
duces the decoherence associated to many photon entan- 
glement. Moreover, an striking consequence of such a 
higher dimensional encoding is that violation of local re- 
alism for two maximally entangled qudits is stronger than 
for two maximally entangled qubits, and increases with 
d 0,^3- For quantum computation applications, OAM 
photon eigenstates could even enable to optimize certain 
quantum computing architectures, where a compromise 
between the number of required qubit and qudit states 
exists |T^. 

The aim of this paper is to develop a general quanti- 
zation scheme of field operators in both the nonparaxial 
and paraxial regimes of light propagation. Within the 
nonparaxial regime, the obtained operators posses the 
suitable phase structure that will straightforwardly allow 



to proceed, at a later stage, to paraxial field operators. 
In this regime, the field operators can conveniently be 
expressed in terms of eigenstate modes of the paraxial 
orbital and spin angular momentum operators. Our ap- 
proach is further extended to the optical phase space. 
The paper is organized as follows: Section II reexamines 
the problem of the separation of the angular momentum 
of a classical electromagnetic field into orbital and spin 
angular momentum components and provides a general 
approach as to whether such a decomposition is phys- 
ically meaningful. Section III adapts the previous ap- 
proach into the field quantization formalism. Section IV 
presents a powerful scheme, based on the Wigner repre- 
sentation, to describe geometric transformations of pho- 
tons prepared in states bearing OAM. Conclusions of the 
paper are drawn in Section V. 



II. CLASSICAL APPROACH REVISITED 

Energy, linear momentum and angular momentum 
constitute the key physical quantities that characterize 
the electromagnetic field configurations. First of all, they 
are constants of motion. Their conservation can be cast 
as a continuity equation relating a density and a flux 
(tensor) density, or current, associated to the conserved 
quantity. The total free electromagnetic angular momen- 
tum J(ro) in a given volume V with respect to a point i"o 
is defined by [lij 

J(r ) = £ / d 3 r (r - r ) x (E x B) 
Jv 

= J(0)-r o xP, (1) 

where r is the position vector, E and B are the elec- 
tric and magnetic fields, and P denotes the total linear 
momentum of the free electromagnetic field. We shall 
focus our study on J(0) = J. Notice that J is defined 
in an analogous manner as the angular momentum of a 
system of massive particles. It can be expressed as the 
integral of an angular momentum density which is equal 
to the cross product of r with the linear momentum den- 
sity £o(E x B). The conservation of J is guaranteed if the 
flux of angular momentum through the surface S enclos- 
ing V vanishes (e.g. for fields that decay sufficiently fast 
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when V becomes large). This is reflected by the integral 
form of the continuity equation for the i-component of 
the angular momentum 



d 4 = -jM lidSl 



(2) 



where M H = e ljk rj [S M (e E 2 + fi 1 B 2 ) /2 - e E k Ei 
—[-Iq 1 B k Bi\ is the angular momentum flux tensor |63|. 

It is well known that in classical mechanics the total 
angular momentum of a system of point matter particles 
can be separated into two contributions. A first one de- 
scribes the angular momentum associated with the cen- 
ter of mass motion of the whole system, while the second 
one refers to the relative angular momentum of the con- 
stituent particles with respect to the total center of mass. 
The center of mass contribution is thus dependent on the 
choice of the reference frame, whereas the relative part 
is independent. Both quantities obey independent evo- 
lution equations. In quantum mechanics, in addition to 
these two contributions (which one could consider as giv- 
ing rise to a purely orbital angular momentum), there is 
an intrinsic or spin angular momentum (independent of 
the choice of a reference frame) with no classical analog 
for point particles. Now, the situation for the free electro- 
magnetic field is more subtle. A similar separation of J 
into strict orbital and spin angular momentum vectors is 
known to be impossible because no reference frame exist 
for the photon in which it is at rest 0,0,0. It is how- 
ever feasible to decompose J into two observables [l7j . 
whose physical meaning will be provided below, as 



with 



L = e J2 / d 3 r Ej- {r x V) A; 



d 6 r E 



(3) 



(4) 



(5) 



where the symbol _L denotes the transverse component 
of the fields (recall that transverse components of any 
field F satisfy V • Fj_ = 0). Since L and S involve the 
transverse part of the vector potential A, they are gauge 
invariant. Notice that, at variance with L, S is inde- 
pendent of the definition of the origin of the coordinate 
system. As it occurs for J, one may show that both L 
and S satisfy continuity equations similar to Eq. J2J). 

Within a purely classical description of paraxial light 
propagation, the total optical angular momentum along 
the propagation direction can be decomposed into spin 
and OAM contributions, each associated with polar- 
ization and phase distribution of the beam, respec- 
tively 0, (see also the excellent reviews and refer- 
ences therein on this subject [ToL |2C| ). As we will prove 
below, these two contributions correspond to the parax- 
ial versions of L and S along that same direction. Our 
framework is completely general, thus allowing to envis- 
age the proper quantization of the fields at a later stage. 



We begin by considering in the Coulomb gauge the 
well-known expansion in the continuous plane- wave basis 
of the vector potential A (which is henceforth assumed 
to be transverse) 



A(r,i) = J2j 



d 3 k 



(16^ 3 £ c|k|) 



1/2 



[e CT (k) a a (k) < 



( i(k-r— c|k|4) 



C. C. 



(6) 



where ao-(k) are the complex amplitudes corresponding 
to the two circular polarization unit vectors e CT (k) (a — 
+1 for right-handed and a = — 1 for left-handed). They 
satisfy e CT (k) • e*, (k) = S aa r and k-e CT (k) = 0. Expansion 
© is a solution of the d Alembert wave equation. 

Our first aim is to derive an equivalent expression to 
Eq. ijfJJl which could then lead in a natural way to its cor- 
responding paraxial limit. The essential feature of any 
paraxial field is that it can be represented as an envelope 
field modulating a carrier plane wave with wave vector 
k . Without loss of generality we take k = k u z , with 
ko > and unitary vector u z , which implies that the car- 
rier plane wave propagates along the positive z direction. 

Let us introduce the trivial identity 



iko(z-ct) 

dko ,^_^ *[fco-/(k)] = l, (7) 



where the function /(k) > 0, which at this stage can 
be any arbitrary function, will be specified below. Upon 
multiplication of Eq. © by J7J) and rearranging, we have 



/>oo 

A(r,t)= / dfc e lfeo(z - ct) A fco (r,i) 
Jo 



(8) 



Now, by imposing Afc (r, i) to obey the paraxial wave 
equation, we obtain the explicit dependence of /(k) 



/(k) 



k z + ^Jkl + 2<z 2 



(9) 



where we have denoted k = q + k z u z , q and k z be- 
ing the transverse (in the x-y plane) and longitudinal 
(along the z-axis) wave vectors, respectively. Equa- 
tion yields a dispersion relation between ko and k 
which gives rise to the constraint imposed by the ar- 
gument /(k) in the delta function of Eq. J2J). Namely, 

5[k -f(k)] = s[k g -(ko-^)] (1 + tf 2 ), with = 

qj \j2k 2 a parameter that governs the degree of paraxi- 
ality. 

The above transformations enable us to map the vari- 
able k Zl defined in the whole real axis, into the more 
convenient positive variable fco- I n this way, one may 
easily carry the integration with respect to k z and cast 
Eq. © in a form that displays both the paraxial and 
nonparaxial contributions in a more transparent fashion 
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AM) = Y,J o °° dk ° J ^ 



x e lfco(z - ct) exp 



1 + tf 2 



1/2 



167T 3 £ CW1 + ^ 4 

iq • r ± - ik d 2 z - ick Q U/l + tf 4 - 1 j f 
I 



[e a [q,fc (l-tf 2 )] a ff [q,k (l-V 2 )] 
c. c.| . 



(10) 



Isotropy of free space allows us to choose the unit polar- 
ization vectors e a [q, fco(l — z? 2 )] in a variety of ways. In 
what follows, it will be convenient to write them as 



e CT [q,fc (l-i? 2 )] = 



V2 



u, 



(l — i? 2 ) 



2d 2 



1 + t9 4 



(11) 



where u p = q/g, u v = u 2 x q/g, and is the polar an- 
gle in cylindrical coordinates. The circular polarization 
vectors e a [q, fco(l — i? 2 )] are both of them orthogonal 
to q + u z fco(l — 1? 2 ), as they should. It is important to 
emphasize that the field (|10[l is exactly equal to the start- 
ing expansion (JSJ. It still obeys the d'Alembert equation 
for all times t > 0. However, it now exhibits the suit- 
able structure to perform the paraxial approximation for 
which !?< 1. In this limit, Eq. llUfl reduces to 



A P (r,t) = J2 



dkn 



d 2 q 



o (16^ 3 e cfco) 1/2 

x e ik a (z-ct) e i(q-r ± -k -9 2 z) , c c 



e CT a a (q, fc ) 
(12) 



where we have only retained the quadratic dependence 
on i3 in the second phase factor of Eq. (jTU)l as the rel- 
evant paraxial contribution. The polarization vectors 
e a = (u x — iaUy)/y/2 are now independent of q and 
ko (u x and Uj, are the unit vectors along x and y di- 
rections). Notice that the structure of the q-integrand 
on the right-hand-side of Eq. Q12 I resembles the well 
known paraxial angular spectrum [2l| . One may exploit 
this fact and expand, rather than in transverse plane 
wave components, into a new transverse state basis: that 
of Laguerre-Gaussian (LG) modes, LGi yP (r±, z; k ), and 
their Fourier-transformed profiles, ££7z, p (q), at z = (see 
Appendix for details). Using their closure relation, one 
has 

Jq-r i -fe i5 2 zl 



,p(r_L, z; ko) 
i, P 



(13) 



Since the LG modes constitute a complete, infinite- 
dimensional basis for the solutions of the paraxial wave 
equation, any spatial beam satisfying the paraxial wave 
equation can therefore be represented in the LG basis in 
terms of an infinite expansion with complex amplitudes 
av,;,p defined by 



(14) 



cv.z.p (ha) = d q£Gt „(q) a a (q, k ) 



Hence, one may cast Eq. (|12[) as 

dko 



A P (r,t) = J2 



o (16^ 3 e cfco) 1/2 



• a<r,z,p (ko) 



<7,l,P ' 

x e iko( - z -^LGi tP (rx,z;ko)+ c. c.l. (15) 

The paraxial electric and magnetic fields follow from 
Eq. Ijl5|l via the well-known relations Ep = —dAp/dt 
and Bp = VxAp. It is now possible to show, by em- 
ploying our convenient representation of these paraxial 
fields into expressions @ and (J5J, that the z-components 
of L and S are given by 

/>oo 

L z = Vl dko \a a , ltP (k ) | 2 , (16) 

poo 

S z = ex / dk |a CT ,z,p (ko) | 2 . (17) 
, Jo 

(7,1, p 

That is, within the paraxial approximation, the total 
angular momentum J z along the beam propagation di- 
rection, i.e. along z, can be decomposed into the so- 
called orbital L z and spin S z angular momenta which 
are related to the azimuthal phase dependence of the 
LG mode basis and their corresponding circular polar- 
ization state, respectively. Equations i|16|) and (|17fl gen- 
eralize the well known results of Allen and coworkers 
with the remarkable feature that they now posses the 
appropriate form to carry out the field quantization. 
Notice also that, since the total energy of any parax- 
ial spatial beam is Hp = Eo J d 3 r [Ep + c 2 Bp] /2 = 
S<r,z,p Jo°° *° c ko\acr,i,p (ko) | 2 , in units of K, the ratio 
L z /Hp can be conceived from a semiclassical point of 
view as the OAM per photon. 



III. PARAXIAL QUANTIZATION 

In the previous section we have derived general classi- 
cal paraxial expressions for the orbital and spin angular 
momenta starting from the continuous plane wave ex- 
pansions of the fields in free space. In this section we 
shall undertake the paraxial quantization of the fields. 
This problem has been considered by several authors in 
the past (see Ref. and references therein). There, 
the approach was approximated, often requiring per- 
turbation expansions which did not provide any clear 
and manageable expressions for the paraxial quantum 
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modes, and, more important, they resorted to the quasi- 
monochromatic approximation which is unsuitable for 
describing such nonclassical processes as the propagation 
of broad-band entangled photons generated by sponta- 
neous parametric down conversion. 

The framework developed in the previous section can 
easily be adapted for field operators. The essential step is 
to substitute the set of complex amplitudes a by creation 
and annihilation operators. At a suitable stage we will 
find more convenient to introduce a different set of cre- 
ation and annihilation operators: those associated with 
the LG modes. 

The expansion for the vector potential operator A(r, t) 
in terms of continuous plane-wave mode operators [recall 



J 



Eq. ®] is 



A(r,t) = E/ d3k 



1/2 



167r 3 £oc|k| 
[e CT (k)a CT (k)e i ( k - r - c l k l t >+ h. c. 



(18) 



The annihilation and creation operators a a and aj. satisfy 

the usual canonical commutation rules [a CT (k), a£.,(k')] = 

<W<$ (3) (k-k'). 

Using again the trivial identity J7J, we can derive the 
quantized version of Eq. (|10|l 



x e ifco(z - ct) exp 



hi + i? 2 



1/2 



167r 3 e cfc \/l + tf 4 
iq • r x - ik i9 2 z - ick (y/l + ?9 4 - 1 j t 

I 



[e a [q,A;o(l-l? 2 )] a a [q, fc (l - tf 2 )] 
h. c.j . 



(19) 



The unit polarization vectors e a [q, feo(l — $ 2 )] are given 
by Eq. 1)11(1. Equation l|19|) . together with its classi- 
cal counterpart l]lUp. constitute the first main results 
of this work. We stress that the field operator (|19fl is 
exactly equal to expansion JTSJ, and thus, it obeys the 
d'Alembert wave equation for any time t > 0. 

We may now build the paraxial version of Eq. I|19fl . As 
mentioned previously, this corresponds to the limit with 
1? <C 1. Recalling the transformation (|13|) for the LG 
modes basis, Eq. (|19|) reduces to the form 



A P (r,t) = 



dk ° I i« • 
\lo7r J £o c « ; o 



1/2 



e CT a<r,z, P (fc ) 

x e 4fco ^" ct )LG; !P (r ± ,z;fc ) + h. c.l, (20) 



where the circularly-polarized polarized vectors e a = 
(u T — i<7\iy)/\f2 are again independent of q and fco- At 
this point, we have introduced the LG mode annihilation 
operators 



a<T,z, P (fco) = y d 2 q££ z * P (q) «<7 (q 5 fc o) , 



(21) 



which, by employing the commutation relations 

[fi<r(q,fo),a£/(q',fco)] = <W<5 (2) (q- q')<K fc o - fco)> and 
the orthonormalization conditions for the LG modes, sat- 
isfy the commutators 



a a ,i, p (ho) (^o) = <W£«'<W<5(feo — k' Q ) . (22) 



We have now at our disposal all the necessary ingredi- 
ents to examine the quantization of angular momentum. 



For massless particles, such as the photon, the only phys- 
ically meaningful component of the total angular momen- 
tum operator J is the one along their propagation direc- 
tion hj. Of course, one may formally quantize the de- 
composition of J = L+S where the operators L and S are 
given by Eqs. Q and © with the classical fields replaced 
by field operators. They can be cast in the form [23| 

1 = -f E /rf 3 k{e* a (k)at(k)e ic l k l*[kxV k ] 
x e^(k)a CT ,(k)e- ic l k l*- h. c.} , (23) 



and 



(24) 



Notice, however, that neither component of L and S is 
a true angular momentum operator since they do not 
fulfill the usual SU(2) commutation relations. Instead, 
they read as: [Si,Sj] = 0, [Li,Sj] = iheijkSk, and 
[Li,Lj] = iHsijk(Lk — Sk)- Nevertheless, it is still pos- 
sible to choose any two commuting components L,; and 
Si, and thus construct simultaneous eigenstates of these 
operators. 

Within the paraxial approximation discussed above, 
we may concentrate ourselves on the two commuting 
operators L z and S z . Upon inverting Eq. I|21l) . yield- 
ing a a (q,k ) = J2i ip £Qi >p (q) a CT) j lP (fc ), we find from 
Eqs. and pijl that the paraxial OAM and spin op- 
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erators are finally given by 

/•OO 

L z = f^y^l I dko a/ al (fco) & ff( j jP (fc ) , (25) 

er,£,p 

/>OC 

5*2 = ^X " / dk oal lp {k )a a: i :P (k ) . (26) 

This means that the most general paraxial one-photon 
state can then be described as consisting of arbitrary su- 
perpositions of eigenstates of L z and S z 



E 



dkoC^M&ltMfl), (27) 



where |0) is the vacuum state. The complex coef- 
ficients C ar i iP (ko) satisfy the normalization condition 
So- 1 p Io° dko\C a j tP (ko)\ 2 — 1, and can be interpreted as 
the probability amplitudes for finding the photon in an 
eigenstate \a,l,p,ko) (Fock state) with circular polariza- 
tion <r, wave vector fco along the z-axis and corresponding 
to a LG mode having indices I and p, that is, with a well 
defined spin and OAM in the direction of the z-axis. 

Interestingly enough, if one uses Eq. (|21|l and rewrites 
the one-photon state l|27|l as 

/■OO P 

\i>) =J2J dk J d 2 q/ <T , Lp (q,fc )4(q^o)|0), (28) 

where / ff ,/, p (q, k ) = Ccr,z, p (fco)£!?i,p(q), one can conceive 
/cr,;,p(q, kg) as the components of the paraxial photon 
wave function in momentum representation. This implies 
that, within the paraxial approximation, the energy den- 
sity of a single photon can be localized in the transverse 
plane orthogonal to their main propagation direction {z- 
axis) with a Gaussian-dependence falloff. Such a spatial 
localization is not in contradiction with the exponential 
(but less than Gaussian) falloff localization limit shown 
by Bialynicki-Birula [24j which applies to three (and thus 
to nonparaxial photons) rather than to two spatial di- 
mensions (see also Refs. [2^, |2f| 153, for closely re- 
lated discussion on this). In fact, the Paley- Wiener the- 
orem still affects the maximum localization of the energy 
density along the propagation direction. Moreover, the 
characteristic length that controls the transverse spatial 
extension of the LG mode cannot certainly take arbitrar- 
ily small values, but those which are compatible with the 
paraxial approximation, i.e. much larger than the photon 
wavelength A = 27r /ko. 

It is important to point out that although the described 
paraxial quantization has been carried out for LG modes, 
which constitute the natural basis for operator L z , it is 
straightforward to show that one could have chosen as 
well other paraxial modes such as the Hermite-Gaussian 
(in Cartesian coordinates) and the Ince-Gaussian |29j (in 
elliptical coordinates). Moreover, our nonparaxial ex- 
pression (|T§|> could readily accommodate photon states 
in Bessel (diffraction-free) modes j23, Ej and oblate 
spheroidal mode solutions of the Helmholtz equation , 
to name just a few examples. 



IV. PHASE SPACE PICTURE OF OAM 
PHOTON STATES 



Phase space, which is a fundamental concept in classi- 
cal mechanics, remains useful when passing to quantum 
mechanics. In a similar fashion with probability density 
distribution functions in classical systems governed by 
Liouville dynamics, quasiprobability distributions have 
been introduced in quantum mechanics. They can pro- 
vide a description of quantum systems at the level of 
density operators (although not at the level of state vec- 
tors). Among them, the Wigner function stands out be- 
cause it is real, nonsingular, yields correct quantum me- 
chanical operator averages in terms of phase space in- 
tegrals and possess positive-definite marginal distribu- 
tions |33l l34l Is^jl . It is, however, only positive for Gaus- 
sian pure states, according to the Hudson-Piquet theo- 
rem. 

By exploiting the analogy between classical and quan- 
tum mechanics with geometrical and wave optics, Wigner 
distributions have been developed in the context of clas- 
sical wave optics of both coherent and partially coherent 
light fields |3y, |33, |38j , where they are Fourier-related to 
the cross-spectral densities. Particularly outstanding has 
been the symplectic invariant approach by Simon and 
Mukunda 39], which has been applied to anisotropic 
Gaussian Schell-model beams via the relation between 
ray-transfer matrices of first-order optical systems and 
unitary (metaplectic) operators acting on wave ampli- 
tudes and cross-spectral densities. 

It is well known that a convenient way to visualize 
the transformation of qubits is provided by the Poincare 
(also known as the Bloch) sphere representation. For po- 
larization states, the north and south poles correspond 
to right- | cr = +1) and left-handed \<r = —1) circularly 
polarized eigenstates, respectively. More generally, any 
completely polarized state can be described as a linear 
superposition of left- and right-handed circular polariza- 
tion in the form (up to a global phase) 



\0,<P) 



cos - \a = +1) 



(29) 



which, on the Poincare sphere, corresponds to a point on 
the surface having angular coordinates 9 and ip. 

In an analogous manner with the above picture for 
polarization states, a Poincare sphere was introduced by 
Padgett and Courtial to represent par axial first-order- 
mode spatial beams carrying OAM [40|. Its underlying 
SU(2) symmetry was subsequently shown |4l|. In this 
picture, the poles of the sphere correspond to LG modes 
with radial-node number p — and topological charge 
I = ±1 (plus and minus standing for the north and south 
poles, respectively). Hence, in complete analogy with 
Eq. I|29|) , any state on the first-order-mode sphere can be 



written as 



>¥>) 



1=1, p=0 



cos- \l = l,p 
e iv sin| |Z = - 



= 0) 

-l,p=0). (30) 



This superposition generates other structurally stable 
states. For example, in the equatorial plane of the sphere, 
combinations with 9 = tt/2 and ip = (9 = n/2, ip = it) 



, n 



) having in- 



give rise to Hermite-Gaussian modes |n a 
dexes n x = 1, n y = (n x = 0, n y = 1). 

Though the first-order orbital Poincare sphere consti- 
tutes an elegant framework to represent families of states 
bearing OAM and their transformation, as points and 
paths connecting these points on the sphere, higher-order 
modes cannot be described by Eq. (J30J. That is, states on 
higher-order Poincare spheres involve more complex su- 
perpositions of LG mode states. A generalization would 
provide, among the many possibilities, direct means to 
visualize the development of geometric phases in optical 
beams |42ll43| and photon states. 

Recently, we were able to carry out the abovemen- 
tioned generalization to all higher-order orbital Poincare 
spheres, and extend the concept of geometric phases in 
paraxial optical beams under continuous mode transfor- 
mations [44|- Via an SU(2) Lie- grou p operator algebra, 
we mapped spin coherent states |45| into families of spa- 
tial modes carrying OAM and belonging to such gener- 
alized Poincare spheres. Remarkably, these families of 
spatial modes could be represented in a compact form 
by resorting to the Wigner function formalism, allowing 
to reveal their hidden symmetries. Our aim here is to 
explicitly show this generalization by constructing sets 
of OAM photon states \8, (p)i tP as points (6, ip) on orbital 
Poincare spheres Oi :P (which can be labeled by I and p or 
equivalently by I and the sphere-order N = 2p+ \ l\ > 0). 
The continuum of states on Oi >p can be generated from 
LG mode states [64[ \l,p) = a\ p \Q) (fixed on the poles of 
the sphere) through the following unitary operations: 

\9,<p) l>p = eM-^t-u v )\l,p) = U(9,<p)\l,p). (31) 

Here, C is a true angular momentum operator (not to be 
confused with L, see below) and u v = (— simp, cos ip, 0) 
a unit vector in the equatorial plane of the sphere. 
The action of the unitary (metaplectic) operator U on 
states \l,p) can be interpreted as a counterclockwise 
rotation of 9 about u v that takes the axis contain- 
ing the poles into the direction with unit vector u r = 
(cos if sin 9, sin if sin 9, cos 9). The important observation 
is that such a rotation gives rise to multidimensional su- 
perpositions of LG mode states that can be cast as 

i P 

\6><p)i, P = E ^2Q',p'(o,v;1,p)\i',p') , (32) 

l' = -lp'=Q 

where 2p' + \l'\ = 2p + \l\ = N. The complex coeffi- 
cients Ci> iP '(6,ip;l,p) depend on the point (9,<p) on the 




FIG. 1: Orbital Poincare sphere of second-order-modes. The 
poles of the sphere correspond to Laguerre-Gaussian modes 
with I = ±2 and p — (positive and negative signs for north 
and south poles, respectively). The states in the equatorial 
plane (9 = it/2) with tp = 0, and ip = n yield the Hermite- 
Gaussian modes having indexes n x = 2, n y = 0, and n x — 0, 
n y = 2 (not shown), respectively. 



iVth-order Poincare sphere 0^ p . One may easily ver- 
ify that Eq. 13211 includes the particular case of Eq. H3U|) 
with Ci,o(0,p;l,O) = cos 6>/2, C o ,o(0,p;l,O) = and 
C-ifi(9, tp\ 1,0) = e ,¥, sinfl/2. Figure ^ depicts represen- 
tative modes associated to their angular orientations on 
the second-order-mode Oi :P sphere (N = 1 = 2). 

The (9, ^-distribution of states IpH)) on Oi iP can be 
obtained by employing the Wigner function representa- 
tion. In the optical phase space, let = (x, y) and 
P = (Px,Py) denote the transverse position and mo- 
mentum (normalized wave vector) variables, respectively, 
and fj_, p be the associated canonical Hermitian op- 
erators. The only nonvanishing commutation relations 
among these operators are [ijpj = [y,Py] = iA. The re- 
duced wavelength A = \/(2ir) = 1/ko plays here the op- 
tical analog of h. For convenience, we arrange the phase 
space variables and the canonical operators in column 
vectors C = (x,y,p x ,P v ), and C = (x,y,p x ,p y ). In terms 
of £, the components of the operator C are 



r — 


x 2 -f (p\ 


2wl 


r — 

J^y — 


Xy PxPyUll 


wl 4A 2 ' 


L - 


zpy - yf>x 

2A 



8A^ 



(33a) 
(33b) 
(33c) 



They satisfy the usual SU(2) angular momentum com- 
mutators [Ci,Cj] — itijkCk (at variance with the com- 
ponents of operator L). Of these SU(2) generators, only 
C z represents real spatial rotations on the transverse x-y 
plane, whereas C x and C y represent simultaneous rota- 
tions in the four-dimensional phase-space: C x produces 
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rotations in the x-p x and y-p y planes by equal and op- 
posite amounts, whereas C y gives rise to rotations in the 
x-py and y-Px planes by equal amounts. Hence, only C z 
does actually correspond (it is proportional) to compo- 
nent L z of Eq. Ij25(l . There is, however, no analog corre- 
spondence between C x and C y and the x and y compo- 
nents of the operator L, respectively. Operators C x and 
C y necessarily involve a change of the spatial modes. 

The Wigner representation of a photon in a pure state 
V'(P) = (Pl^) is 

1 f°° 



x ^(p + ~OV>*(p-~£)- 



(34) 



As shown recently in Ref.|ij, it is possible to obtain the 
Wigner representation of states \9,(p)i iP without explic- 
itly calculating the integrals in Eq. 1)340. To this end, 
we invoke two remarkable properties: (i) On account 
of the Stone-von Neumann theorem, unitary operators 
whose generators are quadratic in £ [such as U(9, ip)] 

induce linear canonical transformations, T : £ — * T£, 
in the optical phase space; (ii) under the action of T 
the Wigner function experiences a simple point trans- 
formation W(£) -> W'(C) = WiT^Q In our 
case, the linear canonical transformation generated by 
the quadratic operators l)33|) results from the relation 
TC = U^GU, and reads as 



T = 



( 


co 










ce 




se s v 


se c v 


V 


zo 

se c v 


se^ v 


zo 


zo 



-z se s v z se c v \ 

z s g c v z se s v 
c e 
c e J 



(35) 



where ce — cos(#/2), sg — sin(6*/2), c v — cosip, s v = 
sirup (recall that zq — w 2 /(2X) is the Rayleigh range and 
wq the mode width at z = 0). Notice that T has the 
form of a symplectic ray-transfer matrix of a generally 
anisotropic first-order system [46l |47| ; it is an element of 
the symplectic group Sp(4, R), that is, det T = 1, and, 
under transposition, TAT* = T'AT = A, where A is a 
real antisymmetric nonsingular 4-dimensional symplectic 
metric matrix 



A 



02x2 
~ 12x2 



12x2 
02x2 



(36) 



Also, the action of T is independent of the chosen states 
at = ip = 0, that is, besides \l,p), one could have chosen 
other state vectors, for instance, Hermite-Gaussian states 
\n x , n y ) with sphere-order N = n x + n y = 2p + 

The key point is thus to observe that owing to the 
unitary relation (|31|l between states belonging to the 
same sphere Oi jP , knowledge of the Wigner function of 
any given state on Oi jP allows one to determine the 
Wigner function of all states on that same sphere. LG 
states constitute the convenient choice here. Using their 



Wigner representation l48j . together with property (ii) 
and Eq. I|35() . the found normalized Wigner function 
is 



(-1) 



iv 



7T2A 2 

x Ln+i (Qq + 4Q • u r ) 



Qo L»^ (Q -4Q-u r ) 



(37) 



where Q = 2[x 2 + y 2 + {p\ +p 2 y )z 2 ]/w 2 , L m {rj) are the 
mth order Laguerre polynomials, and the quadratic poly- 
nomials Q(C) = {Qx, Q y , Qz) follow from C x , C y , and C z 
in Eqs. Q by replacing C -> C- When 9 = (6 = n) 
one recovers from Eq. I)37|l the Wigner function of LG 
states \l,p) (| - l,p)). If 9 = tt/2 and ip = (9 = tt/2 
and <p = 7r) one obtains the Wigner function of Hermite- 
Gaussian states | rise, N — n x ) [\N — n v , n y )). 

Equation 1)37(1 is a strictly positive and angle- 
independent Gaussian distribution only when I = p = 
(in this case its associated Poincare sphere becomes de- 
generated, i.e. all points (9, tp) on the sphere repre- 
sent the same Gaussian mode state). Moreover, though 
Wi, p (G;9,ip) does not explicitly contain the propagation 
variable z, its spatial evolution along z can be fully de- 
scribed by applying a Galilean boost rj_ — > rj_ — zp. 

The orthogonality relations (scalar product) satisfied 
by states \9, <p)i tP and \9', f')i', p ' are given by the overlap 
integral of their associated Wigner functions 

\l<,p>{8',<p'\0,ip)i,p\ 2 



OO 



J —OO 

\N-l)/2 



(2nXy I d*CWi,, p '(C,6,v')Wi, p (C;9,<p) 

2 



N+l \ / N-l 

k=0 



1 - r 



(38) 



where parameter r = cos 2 [(9 — 9')/2] cos 2 \{<p — >p')/2] + 
cos 2 [{8 + 9')/2] sin 2 [(<p - tp')/2] (notice that < r < 1). 
Equation l|38|l implies the following: (i) any two states 
belonging to different spheres are mutually orthogonal; 
(ii) if I > and p = 0, only states corresponding to 
antipodal points are mutually orthogonal. However, if 
/ > and p > 0, additional points exist on the sphere 
(apart from the antipodal) where their associated states 
are also orthogonal [e.g. if p — 1, Eq. (|38|l vanishes when 
r = (I + 1)/(Z + 2)]; (iii) when / = 0, antipodal points no 
longer correspond to orthogonal states but to identical 
states. 

The expectation value i tP (9, ip\£.\9, ip)i tP may be eas- 
ily evaluated with the help of the Wigner function 13711 . 
Since the operators (|33|l do not involve products of non- 
commuting canonically conjugated operators, their corre- 
sponding phase space representation in the Wigner- Weyl 
ordering is simply given by replacing £ — > £ in Eqs. 133|l . 
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Therefore, 



(£) = l<p (d,<p\C\d,<p) llP 



I 



(39) 



Via the Heisenberg-Robertson uncertainty relation and 
using Eq. i|39|) we obtain that the variances ACi of oper- 
ators (|33|) satisfy the following inequalities 

1, 



Ad ACj > -\e ijk {C k )\ 



(40) 



In particular, states on the sphere equator (with 9 = tt/2) 
yield AC X AC y > 0. Moreover, the OAM carried by 
states \6,(f)i_ p follows immediately from Eq. The 
key observation is to notice that the operator L z , given 
by Eq. (|25(l . for a given polarization a and wave vector 
fc , is related to operator C z by L z — 2hC z . The result, 
i tP (9,(p\L z \9 1 ip}i tP = IhcosO, has a very simple geometri- 
cal interpretation. It is the projection of the unit vector 
u> corresponding to \9, <p)i, p along the vertical axis of 
Oi, p . Hence, arbitrary states on any sphere bear frac- 
tional OAM in units of K. The limiting cases, being rep- 
resented by LG and Hermite-Gaussian states, have the 
well known IH and zero values, respectively [llll9|. 



V. CONCLUSIONS 

We have presented a detailed description of classical 
optical beams and single photon states bearing OAM. 
Our quantum field theory formalism generalizes previous 
studies on this subject and, via phase space methods, 
highlights the inherent symmetries of unitarily related 
families of paraxial spatial modes, of which, those car- 
rying integer OAM (in units of K) constitute one par- 
ticular subset. The A-th order orbital Poincare sphere 
representation enables to visualize mode transformations. 
These transformations correspond to first-order optical 
systems with symplectic ray-transfer matrix T [given by 
Eq. i|35|l ] and evidence that it is possible to manipulate 
single photons prepared in superpositions of OAM states, 
and, thus, implement single qubit gates (which corre- 
spond to particular rotations on the Poincare sphere). 
Combination of these gates with two-photon quantum 
gates exploiting OAM entangled stat es g enerated from 
parametric down conversion [a, ml r7L la. Hoi l49l l50l l5ll l5^| 
constitutes a realistic and a fascinating possibility to- 
wards quantum computation with linear optical net- 
works |53,|5J,|53,|5a,|53. One could also exploit the mul- 
tidimensional Hilbert structure of OAM photon states for 
quantum key distribution (QKD) protocols. At variance 
with QKD schemes employing polarization states, which 
only allow transmission of one key bit per photon and 
require the reference frames of the sender and receiver 
to be aligned with each other [5^, OAM photon states 
are invariant under rotations along their propagation di- 
rection. The continuous reference frame alignment moni- 
torization for polarization states may not seem to strong 



of a restriction for ground-based stations, but it could 
be an important limitation on a moving station such as 
a satellite. The distortions created by atmospheric tur- 
bulence on OAM photon states could be corrected using 
two-dimensional filtering techniques that have been pro- 
posed for image recovering (adaptive optics) under vari- 
ous degradation mechanisms [59j. To conclude, we also 
wish to point out that detection of phenomena related to 
highly energetic photons |60| (X-rays and Gamma-rays) 
carrying OAM could be of interest for astrophysical mea- 
surements of distant cosmic entities (pulsars, quasars, 
black holes) through Compton scattering experiments. 
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APPENDIX A 

For completeness, we provide here a resume of the main 
expressions and properties of Laguerre-Gaussian (LG) 
modes (see also Ref. lolf) . Starting from Maxwell's equa- 
tion in vacuum, one obtains the scalar wave equation 



V 2 F = 



1 d 2 F 



(Al) 



for any of the field components F(r,t) (potential vec- 
tor, electric and magnetic fields, etc). The paraxial ap- 
proximation assumes that if, under propagation along 
one given direction (e.g. the z-axis), the field compo- 
nents evolve in an essentially plane-wave fashion mod- 
ulated by some slowly- varying-amplitude u(r), so that 
F(r,t) — u(r)e lfc °( z-c *) (ko is the wave vector along the 
z-axis), then, u(r) is a solution of the paraxial wave equa- 
tion 2iko(du/dz) + V^u = 0. Depending on the par- 
ticular geometry considered, one may distinguish several 
families of complete, orthogonal-set of solutions which in- 
clude the well-known Hermite-Gaussian beams (in Carte- 
sian coordinates), the LG beams (in cylindrical coordi- 
nates) and the more recently discovered Ince-Gaussian 
beams |29j (in elliptical coordinates). 

LG modes satisfy the paraxial wave equation, which, 
in cylindrical coordinates, reads as 



2ik 



du 



d 2 u 
dr 2 



1 du 
r dr 



1 d 2 u 

r 2 d(f> 2 



(A2) 



Their whole z-propagating normalized profiles are 
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LGi tP (r, 4>, z;k Q ) 



I 2p\ 1 / V2r 
tt(|/| +p)\ w(z) [ w{z) 



2r 2 



w 2 (z) y 



exp 



w 2 (z) t 



exp 



il<j) ■ 



2i?(z) 



(A3) 



where w(z) — woy/l + (z/zo) 2 with wq being the width 
of the mode at z = 0, = z [l + (z /z) 2 ] is the phase- 
front radius, z = k w 2 /2 the Rayleigh range, <&q(z) = 
— (2p + \l\ + 1) arctan(z/z ) the Gouy phase [63, and 
Lp(x) are the associated Laguerre polynomials 



rn— 



(|i|+p)i 



(p — m)\ (|Z| + to)! to! 



(A4) 



The indices Z = 0, ±1, ±2, . . . and p = 0, 1, 2, . . . cor- 
respond to the winding (or topological charge) and the 
number of nonaxial radial nodes of the mode. The wave- 
front (equal phase surface) forms in space part of a heli- 
coidal surface given by l(p + kz =const. The topological 
charge attributed to this wavefront manifold is positive 
(I > 0) for right-handed helicoids, and vice versa. 



J 



r 



Of equal importance are the normalized Fourier- 
transformed LG modes, which, at z — 0, are given by 



£Gi, v {p,(p) = 



x exp 



2n(\l\+p)\ 

il(f 



1 2 



w 2 p 2 



exp 



*^P+\l\) 



(A5) 



where p and tp denote the frequency-space cylindrical 
coordinates. They fulfill the orthogonality conditions 
J d 2 <\CQ1 p {q)CGv tP i(q) = Su'6 pp r, With the help of their 

closure relation £V CGi tP (q)jO,Gi, P (q!) = <5 (2) (q-q') one 
has the following relation 



^£^,p(q)^i,p(q') 



l.„ 



i,p i, P 



where in the last step we employed the propagation for- mula of the angular spectrum pH . This proves Eq. . 
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